The Case for Axion Dark Matter 
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Dark matter axions form a rethermalizing Bose- Einstein condensate. This provides an opportunity 
to distinguish axions from other forms of dark matter on observational grounds. I show that if the 
dark matter is axions, tidal torque theory predicts a specific structure for the phase space distribution 
of the halos of isolated disk galaxies, such as the Milky Way. This phase space structure is precisely 
that of the caustic ring model, for which observational support had been found earlier. The other 
dark matter candidates predict a different phase space structure for galactic halos. 
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One of the outstanding problems in science today is 
the identity of the dark matter of the universe [ij. The 
existence of dark matter is implied by a large number 
of observations, including the dynamics of galaxy clus- 
ters, the rotation curves of individual galaxies, the abun- 
dances of light elements, gravitational lensing, and the 
anisotropics of the cosmic microwave background radia- 
tion. The energy density fraction of the universe in dark 
matter is 23%. The dark matter must be non-baryonic, 
cold and collisionless. Cold means that the primordial 
velocity dispersion of the dark matter particles is suffi- 
ciently small, less than about 10^^ c today, so that it 
may be set equal to zero as far as the formation of large 
scale structure and galactic halos is concerned. Colli- 
sionless means that the dark matter particles have, in 
first approximation, only gravitational interactions. Par- 
ticles with the required properties are referred to as 'cold 
dark matter' (CDM). The leading CDM candidates are 
weakly interacting massive particles (WIMPs) with mass 
in the 100 GeV range, axions with mass in the 10~^ eV 
range, and sterile neutrinos with mass in the keV range. 
To try and tell these candidates apart on the basis of 
observation is a tantalizing quest. 

In this regard, the study of the inner caustics of galac- 
tic halos may provide a useful tool 0, An isolated 
galaxy like our own accretes the dark matter particles 
surrounding it. Cold collisionless particles falling in and 
out of a gravitational potential well necessarily form an 
inner caustic, i.e. a surface of high density, which may be 
thought of as the envelope of the particle trajectories near 
their closest approach to the center. The density diverges 
at caustics in the limit where the velocity dispersion of 
the dark matter particles vanishes. Because the accreted 
dark matter falls in and out of the galactic gravitational 
potential well many times, there is a set of inner caustics. 
In addition, there is a set of outer caustics, one for each 
outflow as it reaches its maximum radius before falling 
back in. We will be concerned here with the catastrophe 
structure and spatial distribution of the inner caustics of 
isolated disk galaxies. 

The catastrophe structure of inner caustics depends 
mainly on the angular momentum distribution of the in- 



falling particles There are two contrasting cases to 
consider. In the first case, the angular momentum dis- 
tribution is characterized by 'net overall rotation'; in the 
second case, by irrotational flow. The archetypical exam- 
ple of net overall rotation is instantaneous rigid rotation 
on the turnaround sphere. The turnaround sphere is de- 
fined as the locus of particles which have zero radial ve- 
locity with respect to the galactic center for the first time, 
their outward Hubble flow having just been arrested by 
the gravitational pull of the galaxy. Net overall rotation 
implies that the velocity field has a curl, V x w 7^ 0. The 
corresponding inner caustic is a closed tube whose cross- 
section is a section of the elliptic umbilic (-D-4) catastro- 
phe 0, It is often referred to as a 'caustic ring', or 
'tricusp ring' in reference to its shape. In the case of irro- 
tational fiow, V X i7 = 0, the inner caustic has a tent-like 
structure quite distinct from a caustic ring. Both types 
of inner caustic are described in detail in ref. Q . 

If a galactic halo has net overall rotation and its time 
evolution is self-similar, the radii of its caustic rings are 
predicted in terms of a single parameter, called jmax- 
Self-similarity means that the entire phase space struc- 
ture of the halo is time independent except for a rescaling 
of all distances by BH): all velocities by R{t)/t and all 
densities by |j42|- t is time since the big bang. For 
definiteness, R{t) will be taken to be the turnaround ra- 
dius at time t. If the initial overdensity around which the 
halo forms has a power law profile 
5M, 1 



where Mi and 5Mi are respectively the mass and excess 
mass within an initial radius r^, then R{t) oc ^3 + 97 [4]. In 
an average sense, e is related to the slope of the evolved 
power spectrum of density perturbations on galaxy scales 
8!] . The observed power spectrum implies that e is in the 
range 0.25 to 0.35 Q. Thejjrediction for the caustic ring 
radu is (n = 1, 2, 3, .. 
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where Wrot is the galactic rotation velocity. Eq. ( (2) is for 
e = 0.3. The a„ have a small e dependence. However, 
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the a„ oc 1/n approximate behavior holds for all e in the 
range 0.25 and 0.35, so that a change in e is equivalent 
to a change in j^ax- (e,Jmax) = (0.30, 0.180) implies 
very nearly the same radii as (e, jmax) — (0.25, 0.185) 
and (0.35, 0.177). 

Observational evidence for caustic rings with the radii 
predicted by Eq. ^ was found in the statistical distribu- 
tion of bumps in a set of 32 extended and well-measured 
galactic rotation curves [9] , the distribution of bumps in 
the rotation curve of the Milky Way , the appearance 
of a triangular feature in the IRAS map of the Milky 
Way in the precise direction tangent to the nearest caus- 
tic ring [13], and the existence of a ring of stars at the 
location of the second in = 2) caustic ring in the Milky 
Way Uj. Each galaxy may have its own value of jmax- 
However, the jmax distribution over the galaxies involved 
in the aforementioned evidence is found to be peaked at 
0.18. There is evidence also for a caustic ring of dark 
matter in a galaxy cluster [l^ . The caustic ring model 
of galactic halos [7| is the phase space structure that fol- 
lows from self-similarity, axial symmetry, and net overall 
rotation. 

Self-similarity requires that the time-dependence of 
the specific angular momentum distribution on the 
turnaround sphere be given by j^, ^ 
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where h is the unit vector pointing to a position on the 
turnaround sphere, and j{h) is a dimensionless time- 
independent angular momentum distribution. In case of 
instantaneous rigid rotation, which is the simplest form 
of net overall rotation. 



j{n) = jmax nx (z xfi) 



(4) 



where z is the axis of rotation and jmax is the param- 
eter that appears in Eq. The angular velocity is 
UJ — ^^^Y^z. Each property of the assumed angular mo- 
mentum distribution maps onto an observable property 
of the inner caustics: net overall rotation causes the inner 
caustics to be rings, the value of jmax determines their 
overall size, and the time dependence given in Eq. ^ 
causes a„ oc 1/n. 

The angular momentum distribution assumed by the 
caustic ring halo model may seem implausible because 
it is highly organized in both time and space. Galactic 
halo formation is commonly thought to be a far more 
chaotic process. However, since the model is motivated 
by observation, it is appropriate to ask whether it is con- 
sistent with the expected behaviour of some or any of 
the dark matter candidates. In addressing this question 
we make the usual assumption, commonly referred to as 
'tidal torque theory', that the angular momentum of a 
galaxy is due to the tidal torque applied to it by nearby 
protogalaxies early on when density perturbations are 



still small and protogalaxies close to one another 13|, [14 1 . 
We divide the question in three parts: 1. is the value of 
jmax consistent with the magnitude of angular momen- 
tum expected from tidal torque theory? 2. is it possible 
for tidal torque theory to produce net overall rotation? 
3. does the axis of rotation remain fixed in time, and is 
Eq. ([3]) expected as an outcome of tidal torque theory? 

MAGNITUDE OF ANGULAR MOMENTUM 



The amount of angular momentum acquired by a 
galaxy through tidal torquing can be reliably estimated 
by numerical simulation because it does not depend on 
any small feature of the initial mass configuration, so that 
the resolution of present simulations is not an issue in this 
case. The dimensionless angular momentum parameter 
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where G is Newton's gravitational constant, L is the an- 
gular momentum of the galaxy, M its mass and E its net 
mechanical (kinetic plus gravitational potential) energy, 
was found to have median value 0.05 [15]. In the caustic 
ring model the magnitude of angular momentum is given 
by jmax- As mentioned, the evidence for caustic rings im- 
plies that the jmax-distribution is peaked at jmax — 0.18. 
Is the value of jmax implied by the evidence for caustic 
rings compatible with the value of A predicted by tidal 
torque theory? 

The relationship between jmax and A may be easily 
derived. Self-similarity implies that the halo mass M{t) 
within the turnaround radius R{t) grows as [^]. Hence 
the total angular momentum grows according to 
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where we assumed, for the sake of definiteness, that the 
infall is isotropic and that j(n) is given by Eq. (U). In- 
tegrating Eq. we find 
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Similarly, the total mechanical energy is 
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R{t) df 
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Here we use the fact that each particle on the turnaround 
sphere has potential energy —GM{t)/R{t) and approxi- 
mately zero kinetic energy. Combining Eqs. ([5]), and 
and using the relation R{tY = ^t^GM(t) [3], we 



find 
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(9) 
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For e = 0.25, 0.30 and 0.35, Eq. © implies A/jmax = 
0.281, 0.283 and 0.284 respectively. Hence there is excel- 
lent agreement between jmax — 0-18 and A ~ 0.05. 

The agreement between jmax and A gives further cre- 
dence to the caustic ring model. Indeed if the evidence for 
caustic rings were incorrectly interpreted, there would be 
no reason for it to produce a value of jmax consistent with 
A. Note that the agreement is excellent only in Concor- 
dance Cosmology. In a flat matter dominated universe, 
the value of jmax implied by the evidence for caustic rings 
is 0.27 0,0]. 

NET OVERALL ROTATION 

Next we ask whether net overall rotation is an expected 
outcome of tidal torquing. The answer is clearly no if the 
dark matter is collisionless. Indeed, the velocity field of 
collisionless dark matter satisfies 



dv 
'dt 



dv 
'dt 



{r,t) + {v{r,t)-y)v{r,t) ^ -V(f){r,t) (10) 



where (j){r, t) is the gravitational potential. The initial 
velocity field is irrotational because the expansion of the 
universe caused all rotational modes to decay away [l6| . 
Furthermore, it is easy to show [3] that if V x v = ini- 
tially, then Eq. (|10|) implies V x i7 = at all later times. 
Since net overall rotation requires V x w 7^ 0, it is incon- 
sistent with collisionless dark matter, such as WIMPs or 
sterile neutrinos. If WIMPs or sterile neutrinos are the 
dark matter, the evidence for caustic rings, including the 
agreement between jmax and A obtained above, is purely 
fortuitous. 

Axions iZ-lOl differ from WIMPs and sterile neu- 
trinos in this respect. Axions are not collisionless be- 
cause they form a rethermalizing Bose-Einstein conden- 
sate [llj. Bose-Einstein condensation (BEC) may be 
briefly described as follows: if identical bosonic particles 
are highly condensed in phase space, if their total num- 
ber is conserved and if they thermalize, most of them 
go to the lowest energy available state. The condensing 
particles do so because, by yielding their energy to the 
remaining non-condensed particles, the total entropy is 
increased. In the case of cold dark matter axions, ther- 
malization occurs because of gravitational interactions 
between the low momentum modes in the axion fluid. 
This process is quantum mechanical in an essential way 
and not described by Eq. (fTO|) . 

Axions form a rethermalizing Bose-Einstein conden- 
sate when the photon temperature reaches of order 100 
eV(//10^^ GeV) [2l[ where / is the axion decay constant. 
By rethermalizing we mean that thermalization rate re- 
mains larger than the Hubble rate so that the axion state 
tracks the lowest energy available state. The compres- 
sional (scalar) modes of the axion field are unstable and 
grow as for ordinary CDM, except on length scales too 



small to be of observational interest [2l| . Unlike ordinary 
CDM, however, the rotational (vector) modes of the ax- 
ion field exchange angular momentum by gravitational 
interaction. Most axions condense into the state of low- 
est energy consistent with the total angular momentum, 
say L — Lz, acquired by tidal torquing at a given time. 
To find this state we may use the WKB approximation 
because the angular momentum quantum numbers are 
very large, of order 10^^ for a typical galaxy. The WKB 
approximation maps the axion wavefunction onto a flow 
of classical particles with the same energy and momen- 
tum densities. It is easy to show that for given total 
angular momentum the lowest energy is achieved when 
the angular motion is rigid rotation. So we find Eq. Q to 
be a prediction of tidal torque theory if the dark matter 
is axions. 

Thermalization by gravitational interactions is only ef- 
fective between modes of very low relative momentum 
2l|. After the axions fall into the gravitational poten- 
tial well of the galaxy, they form multiple streams and 
caustics like ordinary CDM . The momenta of parti- 
cles in different streams are too different from each other 
for thermalization by gravitational interactions to occur 
across streams. The wavefunction of the axions inside 
the turnaround sphere is mapped by the WKB approxi- 
mation onto the flow of classical particles with the same 
initial conditions on that sphere. The phase space struc- 
ture thus formed has caustic rings since the axions reach 
the turnaround sphere with net overall rotation. The ax- 
ion wavefunction vanishes on an array of I lines, which 
may be thought of as the vortices characteristic of a BEC 
with angular momentum. However, the transverse size of 
the axion vortices is of order the inverse momentum asso- 
ciated with the radial motion in the halo, (mw^)^^ ^ 20 
meters for a typical value (10~^ eV) of the axion mass. 
In a BEC without radial motion the size of vortices is of 
order the healing length [i^, which is much larger than 

One might ask whether there is a way in which net 
overall rotation may be obtained other than by BEC of 
the dark matter particles. I could not find any. Gen- 
eral relativistic effects may produce a curl in the velocity 
field but are only of order (w/c)^ ~ 10~^ which is far too 
small for the purposes described here. One may propose 
that the dark matter particles be coUisionfuU in the sense 
of having a sizable cross-section for elastic scattering off 
each other. The particles then share angular momentum 
by particle collisions after they have fallen into the galac- 
tic gravitational potential well. However, the collisions 
fuzz up the phase space structure that we are trying to 
account for. The angular momentum is only fully shared 
among the halo particles after the flows and caustics of 
the model are fully destroyed. Axions appear singled out 
in their ability to produce the net overall rotation implied 
by the evidence for caustic rings of dark matter. 
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SELF-SIMILARITY 

The third question provides a test of the conclusions 
reached so far. If galaxies acquire their angular momen- 
tum by tidal torquing and if the dark matter particles 
are axions in a rethermalizing Bose-Einstein condensate, 
then the time dependence of the specific angular momen- 
tum distribution on the turnaround sphere is predicted. 
Is it consistent with Eq. (j3|)? In particular, is the axis of 
rotation constant in time? 

Consider a comoving sphere of radius S(t) = Sa{t) cen- 
tered on the protogalaxy. a{t) is the cosmological scale 
factor. S is taken to be of order but smaller than half the 
distance to the nearest protogalaxy of comparable size, 
say one third of that distance. The total torque applied 
to the volume V of the sphere is 

f{t)^ d^r 6p{r,t) fx {~V(j>{f,t)) (11) 
Jv{t) 

where Sp{f, t) = p(r, t)^po{t) is the density perturbation. 
Pa(t) is the unperturbed density. In the linear regime of 
evolution of density perturbations, the gravitational po- 
tential does not depend on time when expressed in terms 
of comoving coordinates, i.e. (j){f = a{t)x,t) = 4>{x). 
Moreover 6{f,t) = ^''^^^f^ has the form 5{f = a{t)x,t) = 
a{t)5{x). Hence 

f{t) ^ pn{t)a{tf I (fx 6{x) XX {-V^(f>{x)) . (12) 
Jv 

Eq. ([T2l) shows that the direction of the torque is time 
independent. Hence the rotation axis is time indepen- 
dent, as in the caustic ring model. Furthermore, since 
Po(t) oc a{t)^^, T(t) oc a{t) oc t3 and hence lit) oc L{t) oc 
t3. Since R{t) oc ta + s?, tidal torque theory predicts 
the time dependence of Eq. (jS]) provided e = 0.33. This 
value of e is in the range, 0.25 < e < 0.35, predicted 
by the evolved spectrum of density perturbatuions and 
supported by the evidence for caustic rings. So the time 
dependence of the angular momentum distribution on the 
turnaround sphere is also consistent with the caustic ring 
model. 

In conclusion, if the dark matter is axions, the phase 
space structure of galactic halos predicted by tidal torque 
theory is precisely, and in all respects, that of the caus- 
tic ring model proposed earlier on the basis of observa- 
tions. The other dark matter candidates predict a differ- 
ent phase space structure for galactic halos. Although the 
QCD axion is best motivated, a broader class of axion- 
like particles behaves in the manner described here. 
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